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Abstract: An elementary proof of the 3-body Efimov effect is provided in the case of a 
separable 2-body potential which binds at zero energy a light particle to a heavy one. The 
proof proceeds by two steps, namely i) a, projection of the Hamiltonian in a subspace and 

> . 

, the observation that the projected Hamiltonian generates an arbitrarily large number of 
O ■ 

O I bound states, and ii) a, use of the Hylleraas-Undheim theorem to recover the unprojected 
^ ■ Hamiltonian. The definition of the projectors we use can include mean field distortions. 



O ' 1. Introduction 
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The Efimov effect consists in the claim that those 2-body forces which marginally 
bind pairs of particles in the 2-body problem induce an arbitrarily large number of bound 
states for the 3-body problem. This interesting phenomenon has been studied by many 
authors^"^-*, often in the framework of Faddeev equations, and, in particular, was validated 
by a formal proof^-* . An investigation of the same effect for the 4-body problem concluded^-* 
that the effect does not exist in the case of identical bosons, but is present if three heavy 
particles meet with a fourth lighter particle^'^^. 

In this work, we follow the argument of Fonseca, Redish and Shanley"^'^-* and give 
again a proof of the effect, by means of a very elementary derivation. Symmetrization 
constraints between the three particles are strongly relaxed, two particles being chosen 
as heavy and the third particle as light. We do not use Faddeev equations, but rather 

1 



a Born-Oppenheimer (BO) approach. This is the subject of Section 2. Then Section 3 
reinterprets the BO method as a generator coordinate method, thus as a projection into 
a trial subspace. We complete the proof hj i) estimating the non adiabaticity corrections 
and ii) considering the residual coupling effects of the full Hamiltonian. Self consistent 
distortions of pair wave functions due to the presence of the third particle are discussed in 
Section 4. Finally we provide a discussion and conclusion in Section 5. 

2. Elementary Proof 

We consider the familiar Jacobi coordinates x for pair (12) and y for particle 3 with 
respect to the center of mass of (12). Then we consider the 3-body Hamiltonian 

H = -n^A^/M + w{x)-n^A^/{2ii) + X[v{y-x/2) +v{y + x/2)], /x = m/(l+£/2), (2.1) 

where, for the sake of simplicity, the three particles are spinless, the masses Mi and M2 
of particles 1 and 2 are equal. Mi — M2 — M, and the third mass m is much lighter, 
£ = m/M « 1. 

Both potential functions v and w are assumed to have a short or a finite range, but their 
exact forms, local or even non local, are otherwise irrelevant. For the sake of simplicity, we 
assume that w is too weak to create bound states between the two heavy particles. Hence 
the continuum in the three-body problem could display two sorts of cuts at most, namely 
i)the three-body continuum at energy zero, naturally, and ii)two body-continuum(s) if the 
interaction v were able to bind a (m + M) pair. Actually, the strength parameter A has 
that critical value which makes semi-positive definite the 2-body Hamiltonian 

ho = -h^Az/{2ij) + Xv{z). (2.2) 

Here the degree of freedom z may be either i) z = y or ii) z = y±x/2, but it will be 
understood that /iq acts only upon the degree of freedom y. The fact that /iq generates a 
"zero energy ground state" does not depend on x, which plays for ho the role of a simple 
parameter, just defining the center of the critical potential \v. 
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It is noted that the critical value of A is linked here, from Eq.(2.2), to the value of the 
reduced mass jJL rather than the slightly smaller reduced mass of a pair, mM/ [M + m) = 
m/(l +e). With such a value of A, critical for no bound state is induced for the 2-body 
problem with m/(l + e). The 2-body continuum in the 3-body problem starts at the same 
energy {E = 0, indeed) as the 3-body continuum. At most a resonance may occur in the 
2-body channel at a slightly positive energy E^. 

For the sake of simplicity, the argument which follows assumes that v is a scalar, non 
local, rank one, separable potential^) , 

<p'\v\p>=-f{p')f{p), (2.3) 

where p is the momentum conjugate to z, and p,p' are the lengths of p respectively. 
The well-known binding condition for separable potentials then reads, at the critical value 
of A, 

1/A = I dp[f{p)Y/p^ (2.4) 
Moreover an exactly soluble model is obtained if we choose the following form factor /, 

/(p) = l/[7r(pV7')], (2.5) 

where I/7 defines any suitable short range in coordinate space. With Eq.(2.5), the condi- 
tion, Eq.(2.4), for a critical value of A becomes 

A = 7^ (2.6) 

In the following we set % — 1 and ?i^/{2ii) — 1, for a simpler system of units. We also 
choose p to be the momentum conjugate to y, and introduce the translation operator 
T = exp(zp • x/2) and its inverse = exp(— • x/2). Then that part hf oi H which acts 
upon y reads, 

hf=p'^ + \{TvT-^ + T-^vT) = / - A(T|/ >< /[T"^ + T" V >< f\T). (2.7) 

Parametrized by x, this Hamiltonian hf drives the BO dynamics of the fast, light particle. 
When compared with the semi-positive definite /io, this Hamiltonian hf contains one addi- 
tional attractive potential, and thus induces a truly bound state (square integrable) Xxiv) 
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at a definitely negative energy rjg. Since threshold singularities of a square root nature with 
respect to rj are expected, we define such a square root cohy rjg = — [cog]'^. Naturally o;^ — > 
when X ^ oo. The subscript x will be omitted in the following, unless it is essential. 

The wave function x is easily obtained in momentum representation, according to 

X(p) = A < f\T\x > 2cos(p • x/2) f{p)/{p^ + cj^), (2.8) 

where we take advantage of the symmetry of the ground state, < f\T~'^\x >=< f\T\x > ■ 
Projecting Eq.(2.8) against < /|T, we obtain the equation which solves for the binding 

energy, 

1/A = /(a;,f) +/(a;,0), (2.9) 

with 

I{io,x) ^< /|(/+a;2)-iTV > . (2.10) 
We note incidentally that the condition for the value of A to be critical reads 

1/A = 7(0,0). (2.11) 

With the choice of f{p) as a scalar /(p), the integral /, Eq.(2.10), does not depend 
on the orientation x of x, but only on its length x. More precisely, it becomes 

/•OO 

I{oj,x) = Att/x / dpp sm{px) [/(p)]V(P^ + ^^)- (2-12) 

For the soluble model provided by the choice, Eq.(2.5), a straighforward contour integration 
gives 



/OO 
dp p sin(pa;)/[(p^ + a;^)(p^ + 7^)^] 
-00 



(2.13) 

= [27exp(— a;a;) + {co'^x — ■y'^x — 2-f) exp(— 7a;)]/[7a;(7^ — c<;^)^]. 
It is easy to obtain the value of / for a; = 0, 

/o(a;) = /(a;, 0) = 1/(7(7 + ^)']- (2-14) 

and then, for Eqs.(2.11,6), the number /qo = -^(0,0) = 1/7^. 
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If ujx has a finite fimit c when x — > oo, it is also easy to find that I{u},x) boils down 

to 



= ^'rr-' (2.16) 



2 exp(— Ci;a;) 
a;(7^ — a;^)^ 

This simplification is also useful for the calculation of derivatives of /, because the con- 
tributions of the discarded term, proportional to exp(— 70:), are negligible for deriva- 
tives as well. Let us test this conjectured finite limit of ujx, by means of an ansatz 
uj = c/x + d/x"^ -\- 0{l/x'^) for Eq.(2.9). We obtain the condition 

^ = + ^ + 0{lx^). 2.16 

7^ xYi^ -{c + d/xf/x^f 7[7 + (c + c//a;)/x]' 

The ansatz is then consistent if c and d obey the conditions 

exp(-c) = c, cZ = 3cV[27(l + c)], (2.17) 

hence c ~ 0.5671, d ~ 0.3079/7. We conclude that the "fast" BO Hamiltonian^'^) hf binds 
the light particle at energy 'q{x) = —(? jx^ -\- OiXjx^^ when x — > 00. 

For a better understanding of this l/x^ binding, a thorough examination of Eq.(2.9) 
is in order. We rephrase it as 

l/A-/(a;,0) = /(a;,a;), (2.18) 

and notice on one hand that, because of the Fourier transform introduced explicitly by 
the translation in the definition of /, Eq.(2.10), we can take advantage of Lebesgue's 
theorem to predict a 1/a; factor in the right-hand side of Eq.(2.18) when x — > 00. On the 
other hand, the left-hand side of Eq.(2.18), once it is written as a function of lo rather 
than 77, is a regular function near a; = 0. Since 1/ A = 1(0, 0) , the difference 1(0,0)— /(a;, 0) 
is proportional to first order with respect to oj. Hence the possibility of a connection 
CO oc 1/x. Had A been overcritical, inducing an actual negative eigenvalue —coq, ujq 7^ 0, 
then 1/A = I{ljo,0), and one might have tried a connection uj — lvq oc 1/x, but failed 
because of the resulting decay of exp(— a;oa;) in the right-hand side. This shows the strict 
importance, for the Efimov argument, of the condition that the value A be exactly critical. 



It is now trivial to consider the "slow" BO Hamiltonian, which drives the heavy degree 
of freedom x, 

hs = -/^s/M + w{x) +r]{x). (2.19) 

The effective potential r]{x) induced by the binding of the light particle has the expected 
long range behavior, rj oc —1 / x^ , that justifies the Efimov prediction of an infinite number 
of bound states for the 3-body system. It is interesting to point out, incidentally, that 
the mass ratio s appears strictly nowhere in this derivation of rj. Hence, if the decoupling 
of the full Schroedinger equation into the two separate steps described by hf,hs can be 
justified by any other argument than the Born-Oppenheimer one, it can be claimed that 
the —1/x^ nature of the induced potential is valid also for particles with similar masses. 

We can also offer a qualitative argument for a generalisation of this —l/x"^ nature in 
the case of potentials v other than separable. Let us assume that v is semi-negative definite 
and let p{uj) be the highest eigenvalue of the operator R{uj) = {—vY/'^{jP'-\-uj'^)~^{—vY/'^. 
In the same way, let a{uj, x) be the highest eigenvalue of S{uj, x) = {—v — T^i;T~^)^/^(p^ + 
u)'^)~^ {—V —T'^ vT~'^Y /"^ . The fact that A is critical is expressed by the condition — p(0). 
Then the binding energy in presence of the two potentials TvT~^ and T~^vT is obtained 
when solving for co the "spectral" condition p(0) = a{u!,x). If, when x — > oo, a property of 
the form a{uj, x) = p{uj) + t{ujx) / x occurs, where r is a suitable function depending on the 
product ujx only, then the "spectral" condition which provides uj amounts to p(0) — p{uj) = 
t/x. It is also reasonable to assume that p(0) — p{uj) = 0{u) when u) —>■ 0. Then the whole 
scheme discussed for Eq.(2.18) is recovered, and thus oj ocl/x. 

In summary for this section, elementary arguments indicate that if the distance x 
between two particles is frozen, and if a third particle is "bound at energy zero" to, e.g., 
the first particle, the presence of the second particle will induce true binding of the third 
particle, at an energy rj oc — l/x^ when a: ^ oo. This rj will in turn behave as a long range 
potential, efficient for creating many bound states of the pair made by the first and second 
particles. 
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3. Born-Oppenheimer Approximation and Generator Coordinate Method 



Let (p he a. bound eigenstate of hg, and E the corresponding eigenvalue. For the sake 
of simphcity, we assume in the following that w is a scalar, like rj and v, and we then omit 
vector notations whenever possible. The function </? is a scalar (s state). Although Xx is 
not a scalar, it transforms very simply under any rotation TZ, namely TZ{xx) = Xn{x)- As 
stated already, the vector nature of ^ as a label can often be understood, and/or shortened 
into a simple scalar label x. 

The well-known BO ansatz for the 3-body state ilj{x,y) corresponding to E reads, 

ijj{x,y) ^ (p{x)xx{y)- (3.1) 
This may be written as well under the form, 

i;ix, y)^ j (^(0 5{x - Ox^ (y), (3.2) 
which is nothing but a generator coordinate^^"^^^ expansion on a basis 

(pi{x,y) ^ 5{x - i)xi{y), (3.3) 

of states (j)^ parametrized by the continuous label ^. The "slow" wave function (p{x) is thus 
reinterpreted as a mixture amplitude <^(^) with respect to the generator coordinate ^. 

If this reinterpretation can be shown to be consistent, we will thus obtain that any 
BO eigenvalue E is an eigenvalue of the projection PHP of H in that subspace spanned by 
the wave functions (j)^. This will reinforce the BO result, namely that there are an infinity 
of such negative eigenvalues E. Indeed, we can then take advantage of the Hylleraas- 
Undheim (HU) theorem^), which we summarize by the short statement "the unprojected 
Hamiltonian has more (or at least as many) bound states than (as) its projection into any 
subspace ". 

That is, any exact wave function ^ may be written as 

m = Pm + Q*, (3.4) 
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where Q is a projector orthogonal to P. Then, the effect of the added Q-term modifies the 
spectrum of PHP in such a way that the correct number N of bound states satisfies the 
inequahty 



Here N is the exact number of bound states of the original H, but we have also inserted 
an intermediate number A^^ of bound states in Eq.(3.5) to include the situation where Q is 
not the full complement 1 — P of P, but only part of it, for a variational theory with trial 
states in an enlarged subspace. In many of the earlier treatments of the Efimov effect in 
terms of models and simple ansatz on the wave functions, the effect of the Q-component of 
the wave function was not estimated, which made the proof incomplete. The HU theorem 
for the bound states may be re-stated as a deepening influence of the effective Q-space 
potential PV VP < where is the Green's function in the Q-space. Therefore, all 
the bound states produced by PHP will be uniformly pushed down in energy, and possibly 
additional bound states are created by the Q-component. 

We now proceed by proving that the slow Hamiltonian, Eq.(2.19), is formally equiva- 
lent to a Griffin-Hill- Wheeler (GHW) kernel H{^',^) =< (j)^i\H\(f)^ >, hence a representa- 
tion of the projection PHP of H in the generator subspace spanned by the states (j)^. For 
this, we first notice that such states 0^ are not orthogonal to one another if one measures 
their scalar product by < X^'lxc namely an integration upon y only, but become triv- 
ially orthogonal if an integration upon x is also included. Indeed such states are strictly 
localized in x-space and can be normalized according to < >= ^{C' ~ 0- Nothing 

prevents us from normalizing these states to unity in the j/-space, since they are strictly 
bound eigenstates of the fast Hamiltonian hf, Eq.(2.7), 



Np < Nt < N . 



(3.5) 




(3.6) 



We find incidentally, as a consequence of Eq.(3.6), that 




(3.7) 
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and furthermore, as a consequence of Eq.(3.7), that 

u{0 = j dy xdv) [-^i Xdv)] = j dy [V^ xdv)? > 0- (3-8) 

The GHW kernel thus reads 

< 4>i'\H\4)^ >= < (l)^>\[-A^/M + w{x) + hf{x,y)]\(l)^ > 



(3.9) 



Except for a correction proportional to tt(^), this is nothing but hg in coordinate represen- 
tation. The GHW equation, 

J <cP^^\{H-E)\(P^> v^(O = 0, (3.10) 

and the slow BO equation, 

2e 



A^+w{x)+r]{x)-E)(p{x)=0, (3.11) 



l+s/2 

are thus equivalent if this correction brought by u, a repulsive effective potential, can be 
shown to be negligible. For that purpose we restate Eq.(2.8) in the form 

|x>=(<ClC>)-'/'lC>, (3.12) 

with 

IC >= ip' + ojY\T + T-')\f >, < CIC >= 2 < f\ip' + u;Y\l + T')\f > • (3.13) 
In particular, when a; — > oo, 

d 

a; < CIC > = —^[Ioiuj) + Ioo{i^,x)] 

2 ^ 2exp(— wx) 8ujexp{—ujx) (3-14) 

7(7 + 0;)^ (7^—0;^)^ a;(7^ — a;^)^ 

The right-hand side of Eq.(3.14) has a finite limit ^ = 2[1 -|- exp(— c)]/7'^ when a; — > 00, 
hence < C|C >= 2(1 + c)x/{c^^) + 0{x^) when X — > 00. 



Let us denote by |x' >, |C' > the full gradients of |x >, |C >) respectively, with respect 
to X, including the dependence of co on x. It is easy to find that 

IC' >= f (p' + - T-')\f > -2a;(/ + uj^)-\T + T-')\f > V,- (3.15) 

where we know that u = -(c + 2d/x)x/x'^ + 0{l/x^) = -(c + 2d/x)x/x^ + C(l/a:^). 
Here x is the unit vector defining the direction of x, namely x = x/x. Then we have to 
investigate the behavior of 

lexer 



u x'\x' >= (< CIC >)-^ < C'l (i - ^1^) IC' 



(3.16) 



when X — > oo. 



From the linear asymptotic behaviour, < > oc of the square norm of we 
deduce that the length of its gradient, namely of 2 < (^\(^' >, has a finite limit when 
a; — > oo. Indeed, from the derivative of the right-hand side of Eq.(3.14), we find 

< CIC' >= 



exp{—ujx) ^ 4a;exp(— a;a;) ^ 4exp(— a;x) 



X 



1 



xexp{—u!x) exp(— a;a;) 8exp(— a;a;) 24a;exp(— a;x) 
a;(7^ — a;^)^ a;^(7^— a;^)^ (7^ — 0;^)^ x{'j'^ — uj'^)^ 



+ 



'^^7(7 + '^)^ 0)7(7 + 0;)^ 
a;. 



(3.17) 



the limit of which is (1 + c)£/(c7^) = [exp(c) + l]£/7^. It can be concluded that the 
term — (< > / < CIC >)^ in the right hand side of Eq.(3.16) reinforces the potential 
rj — —(? jx'^ by a contribution — £/[(2 + £)x^\. 

We must now evaluate 

< CIC >= 



1 



< f\v\v^ ^^^Y'\^-T%f > -4uj{c + 2d/x) < f\p-x{p^ +u;'^)-^sm{p-x)\f > fx 



+ Su'ic + 2d/xY < f\{p' + u;')-\l + T')\f > /x^ + Oil/x'). 

The matrix element of p ■ xsm(p- x) can be obtained as 

< f\p- xsm{p- f)(p2 + uj^)-^\f >= Im < f\p- xexp{ip- x){p'^ + o;^)"^!/ >= 



(3.18) 



Xm<f\ 



—i — exp(isp- x) 
ds 



{p'+u;^)-' \f>= 



d \ d ^. , ,^ . 
2^ ^/(-,.x)/2U=x. 



(3.19) 
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Then a somewhat tedious calculation, which cannot be reported here in detail, gives the 
contribution of < C'lC' > to w, 

2e < C'lC > 



(l + £/2)<C|C> V2 + £ 



X 



C7 12 + 36c + 45c2 + - 13c^ - 2c^ o/ 
+ T^TT^Z^,. + ) 



(1 + c)x 6(1 + c)3x2 



0.367 ^ 1-98 



2 + £ y V X a;^ 

(3.20) 



It turns out, as a matter of fact, that < C'lC > has a finite limit 1/(27"^) when x oo. 
Combined with the result already found, < Q\C, >— 2{l+c)x/{c^'^) + 0{x^), this asymptotic 
behavior of < C'lC' > induces a repulsive potential of order 1/x, which seems to contradict 
the leading term —(? jx^ derived from the straight BO approximation. 

We notice however that this repulsive potential £C7/[(2 + £)(1 + c)a;], and also the 
minute corrections brought by u at order 1/a;^, are weighted by £, naturally. Hence the 
Efimov phenomenon is still present for e sufficiently small, because a generic potential of 
the form e/x — l/x"^ generates an arbitrarily large number of bound states. 

More precisely, if we collect all relevant terms contributing to Eqs. (3.9,10), the radial 
equation in s— waves to be studied for binding is, for (p = x(f, 

d? (l + £/2) , , C7 A il + e/2)E\ . ^ ^ „ 

(3.21) 

(1 + £/2)c2 1 12 + 36c + 45c2 + c^ - 13c^ - 2c^ 

Ye 4+ 24(1 + c)3 • ^^-^^^ 

Potentials of order 1/x^ when a; ^ oo do not modify our expected conclusion of a large 

number of bound states and can be simply neglected. We also keep in mind that the 1/x 

and l/x"^ forms contributing to Eq.(3.21) are not valid for small values of x, since the 

potentials they represent are actually finite. 

If now we assume for the sake of simplicity that w has a strictly finite range xq, the 
radial equation Eq.(3.21) can be replaced by 
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with 



The hard core boundary condition (p{xo) = removes all effects of the (finite) potentials at 
short distances. This introduction of a hard core is legitimate, for it can only underestimate 
the number of bound states. 

A rescaling x = 4(1 + c)z/{c^) is then suitable to remove the coefficient of the 1/x 
term. The ratio between x and 2; is a fixed number, which depends only upon c and 7, and 
we define a fixed, rescaled hard core radius, zq — c7a;o/[4(l + c)], accordingly. Finally the 
energy scales according to E' = 4(1 + c)^(2 + £)£^/ (c^7^£), and we find it strictly equivalent 
to study the number of bound states for the following problem, 

72 



/(Pi A \ 

[-— + — -E')<p = 0, z>0, <p{0) = 0. 

\ dz"^ z + zo [z + zoY J 



(3.24) 



The positive coefficient A is of order 1/e and becomes infinite when e ^ 0. It is obvious 
that, as long as z + zq < A/2, then 

* <-.7-^. (3.25) 



z + zq {z + zqY - 2{z + zo)^ 
Therefore the problem, 

(-5^-2(771^-^')^ = ''' "^^^^2, ^(O) = 0(A/2) = O, 



(3.26) 



generates a further lower bound for the number of bound states. From the potential 
present in Eq.(3.26), we finally use the trace formula^^^ for an estimate of a number of 
bound states, 

A/2 A 



Nb= dzz —-^ ^ oc A Log(A). (3.27) 

This number Nb does diverge when e — > 0. Q E D 

To summarize this section, the Born-Oppenheimer method is restated as a special 
case of the generator coordinate method. The latter is a projection of the Schroedinger 
equation onto a trial subspace. The proof that there is an arbitrarily large number of 
bound states in a subspace is extended to the full space of wave functions, by means of 
the Hylleraas-Undheim theorem. In the process of relating BO to GCM, however, a long 
range repulsive correction su is found to perturb the static BO potential 77 which creates 
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the Efimov effect. Small values of the mass ratio, the adiabaticity parameter e, are then 
necessary to validate our proof of an arbitrarily large number of bound states for the 
three-body system. 

4. Introduction of Distortions and Mean Field Approximation 

In order to make less mandatory the restriction of e to small values, we suggest a more 
flexible formulation of the BO method. For this, we complete i?, Eq.(2.1), by a constraint 
upon < X >, via an auxiliary, external harmonic oscillator potential {x — ^'^ with a strong 
spring constant, Lagrange multiplier L. This defines the constrained Hamiltonian 

n = P'^/M + L {x-if +w + p^ + \{TvT-^ + T-^vT). (4.1) 

— * 

It is clear that the ground state |^*|* > oiTi. verifies < $|-|^|$|- >— > ^ when L — > oo. For 
finite values of L, however, the freezing of x is implemented by a wave packet, less stringent 
than a 5— function 5{x — 

For simplicity of notations in the following, vectors will again be replaced by scalars. 
The simplest approximate description of the ground state of H, is obtained by means of a 
Hartree ansatz $^(x, y) — '^i{x) 'E.^{y). It is clear that F is not very different from a sharp 
Gaussian centered at ^, and that, however, it incorporates distortions due to the various 
potentials present into the corresponding Hartree equation, 

[P^/M + L {x-^)^ + w + V-r]'] r = 0, (4.2) 

where V is the mean field potential induced by the convolution of X{TvT~^ +T~^vT) with 
the density [S(y)]^. The corresponding Hartree eigenvalue r]' contains a spurious harmonic 
oscillator contribution, to be discarded if physical interpretations are needed. In turn, S is 
the bound state generated for the third particle according to the second Hartree equation. 
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(4.3) 



where is the mean field potential arising from the convolution of \{TvT~^-\-T~^vT) with 
the density [r(x)]^. In so far as F is strongly localized around ^, there is not much diff'erence 
between and \[v{y — ^/2) + v{y + ^/2)]. Hence there is a strong similarity between this 
second Hartree equation, Eq.(4.3) and the fast BO equation driven by /i/, see Eq.(2.7). 
The same similarity holds for ?7"(^) and the static BO potential r]{^). The former, however, 
contains those corrections arising from the differences between F and a 5— function. It is 
stressed that all these corrections and distorsions represent self consistency between x and 
y, and facilitate binding. 

The generator coordinate ansatz, 

I* >= j F(0 1*^ >, (4.4) 
then leads to the usual GHW equation. 



/ 



(< > -E < $^/|$^ >) F{i) = 0, (4.5) 



where now the overlap kernel A/"^/^ =< ^>^/ > differs from 5{^ — ^'). It is rather similar to 
a Gaussian, as an overlap of the quasi-Gaussian wave packets F^ and F^/ . Its inverse matrix 
square root M~^l'^ will be necessary for the usual GCM deconvolution manipulation. The 
slow wave function will not be F, but rather H^I'^F. 

It may then be convenient to write the Hamiltonian kernel tl^'^ =< > under 

the form 

il^,^ =< $^,| [p'^jM + «; + / +^/^ + \{TvT-^ + T-^vT) -U^ > 

, (4.6) 

=< $^.|(PVM + t(;)|$^ > + < ^^,\{\{TvT-^ ^T-^vT)-U^\^^ > +r/"(0A/'e6 
for a hint that rj" will appear as a dynamical effective potential when the usual decon- 
volution M~^l'^iiM~^l'^ is performed. It will be noticed that, according to the Hartree 
definition ofW, the diagonal matrix element < ^^\\\{TvT~'^ -\-T~^vT')—U.^\^^ > vanishes 
identically. Finally it is known^^^ that the same deconvolution removes the zero-point 
kinetic energy which plagues < $^/|P^/M|$^ > . 

All told, after deconvolution, the present GHW equation, which includes dynamical 
distortions, reads as a generalization of the slow BO equation. The self consistency inserted 
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by the Hartree method, via those distortions included in F and S, is expected to bring 
more binding. It should thus increase the number of bound states and allow larger values 
of e to be compatible with a given, large number of bound states. 

5. Discussion and Conclusion 

In this paper, on one hand, we provide another rigorous proof of an Efimov effect. 
But the proven effect is weaker than the full expected effect. Namely, given the initial 
hypothesis that a pair potential is marginally able to bind, an arbitrarily large number of 
three-body bound states is obtained only if a suitable mass ratio is small enough. This 
may be the result of the BO ansatz, Eq.(3.1), and the particular coordinate system chosen. 

On the other hand, several new results were found. For one, the scheme of our proof 
is based on arguments of moderate technicality only, such as the recognition that the BO 
method is but a special case of a projection of the three-body dynamics into a generator 
coordinate subspace. Hence, we were able to take advantage of the Hylleraas-Undheim 
theorem, which states that there are at least as many bound states in the full space of 
wave functions as there are in a subspace. 

We were also able to relax the BO freezing of the heavy degree of freedom into more 
flexible Hartree calculations under harmonic oscillator constraints. This allows more pre- 
cise descriptions of the wave functions of Efimov states, including mean field distortions. 

Despite the slight technical difficulties we had to face when the effective potential, of 
the form — l/x^, became plagued by corrections of the form e/x,we did not find it necessary 
here to introduce the now well known electron translation factor (ETF) correction of the 
BO wave functions for the calculation of the molecular potential^'^"^^^. An introduction 
of this ETF correction is likely to be in order for future stages of the theory only. We 
conjecture that the e/x perturbation may be eliminated when Eq.(3.1) is improved by 
this ETF asymptotically. The problem of ill-behaved boundary conditions raised by the 
conflict between the representation where the two heavy particles are considered close 
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by and the representation where these heavy particles are taken far apart is not new and 
smooth transitions between such representations have been proposed^^^ We can point out, 
however, that each one of the two competing representations provides a complete basis of 
the Hilbert space of bound states. For our theory of Efimov bound states, the conflict 
between these representations is thus tempered. 

Accordingly, we find it reasonable, while odd at first sight, that the reduced mass 
which we use to define the critical value of the potential strength is a three-body reduced 
mass = 2Mm/{2M + m) rather than a pair reduced mass ji' = Mm/{M + m). This 
choice is indeed imposed mathematically by the "fast" BO Hamiltonian, see Eq.(2.7). But 
it may also receive an intuitive, physical interpretation : that BO bound state Xx{y)-, 
which reaches zero binding when a; — > 00, is even under the exchange of the two heavy 
particles and is defined with respect to a fictitious particle, namely the center of mass of 
these heavy particles. In that sense, the critical condition for marginal binding must refer 
to the relation between the light particle and that center of mass with mass 2M. All told, 
for Efimov states, where long range effects are at work, a pair formed by the light particle 
and one of the heavy ones cannot really be isolated from the other heavy particle. 

Finally, our approach can be extended-*^^^ by the consideration of coupling the three 
possible "channels" defined by the three possible pair partitions. A BO treatment, or the 
constrained Hartree(-Fock) generalization advocated in Section 4, can be undertaken for 
each such channel. There is no doubt that the resulting, projected Faddeev equations will 
be driven by long range potentials similar to those found in the one-channel, mathemati- 
cally rigorous argument detailed in Sections 2 and 3. 
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